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I. INTRODUCTION 


Quark and gluon distribution functions (parton distribution functions (PDFs)) inside the 
hadron play a significant role to study standard model and beyond standard model physics 
at high energy colliders. Consider for example, the Higgs m discovery at ,ar g e hadrorr 
colliders (LHC) { 3 ), 4], The gluon-gluon fusion process (gg —> H) contribute significantly 
to the Higgs production cross section at LHC 5J. Similarly, the gluon fusion and gluon 
fragmentation processes contribute dominantly to the production of the heavy quark, heavy 
quarkouium and jet etc. at high energy colliders. Hence, as the hadron collider’s total 
energy becomes higher and higher, the gluon distribution function inside the hadron plays 
significant role to study standard model and beyond standard model physics at high energy 
colliders. 

Unfortunately, the parton distribution function inside a hadron is not calculated so far 
because it is a non-perturbative quantity in QCD and we do not know the solution of the 
non-perturbative QCD yet. Hence the usual procedure in the high energy phenomenology is 
to determine the parton distribution function from a set of experiments at some momentum 
transfer scale and then determine its evolution to other momentum transfer scale by using 
DGLAP evolution equations jg]. Once the parton distribution function is determined in this 
way then one uses it to predict physical observable at other collider experiments. Hence in 
order for this prescription to work it is necessary that these parton distribution functions 
are universal, i. e., they do not change from one collider experiment to another collider 
experiment. From this point of view it is necessary to use the correct definition of the parton 
distribution function at high energy colliders. If one does not use the correct definition of 
the parton distribution function at high energy colliders then one will obtain divergent cross 
section for physical observable at high energy colliders. This can be seen as follows. 

Consider, for example, the hadron production at high energy colliders. If the factorization 


theorem holds [z-JD] then the formula for the hadron production cross section at high energy 
hadronic colliders is given by 

a(AB^H + X) = j dx i J dx 2 J dz fa/A(x 1 : Q 2 ) f b / B (x 2 ,Q 2 ) d(ab cd) D H/c (z,Q 2 ) 

a,b 


(i) 


where f a /A(x,Q 2 ) is the parton distribution function (PDF), a(ab —$■ cd) is the partonic 
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level scattering cross section and Dh/ c (z,Q 2 ) is the fragmentation function (FF). For non- 
hadronic observable at high energy colliders the fragmentation function Dh/ c (z, Q 2 ) is not 
included in eq. (P). In the above equation a,b,c,d = q,q,g where q, q and g are quark, 
antiquark and gluon. 

Hence if one does not use the correct definition of the parton distribution function or 
fragmentation function in eq. (pQ) then any uncanceled collinear or soft divergences in the 
Feynman diagrams in the partonic level scattering cross section a(ab —> cd ) can not be 
absorbed in the definition of the non-perturbative parton distribution function f a /H(x,Q 2 ) 
or fragmentation function DH/ a (x, Q 2 ) in eq. (CQ). Hence it is necessary to derive the cor¬ 
rect definition of the parton distribution function f a /H(x,Q 2 ) and fragmentation function 
D H /a(x,Q 2 ) from first principles. 

A correct definition of a parton distribution function at high energy colliders must satisfy 
the following three properties: 1) it must be consistent with the number operator inter¬ 
pretation in quantum mechanics, i.e., ^(x)ij>(0) or F^ ua (x)F“ u (0 ) must be proportional to 
quadratic power of the quark or (quantum) gluon field, 2) it must be gauge invariant and 3) 
it must be consistent with factorization theorem in QCD. The present definition of the quark 
distribution function at high energy colliders is consistent with these three properties js]. 
dowever, the present definition of the gluon distribution function at high energy colliders 
8| is not consistent with these three properties as we will see below. 

re present definition of the gluon distribution function at high energy colliders is given 


by 


fg/p( x ) = 


1 


2irxP + 


dy e 


- -ixP+y 


[< P\F +fJ,a (0,y~,0 T ) x Vexp[igT (A)b £ dz~Q +b {Q, z~, 0 T )] x F+ a {0)\P >] (2) 


where 


F^(x) = dMx) - d„Ql(x) + gf abc Q b ,(x)Q c Ax), (3) 

T^ a = —if abc a nd Q fia (x) is the quantum gluon held. 

The expression Vex\)[igT {A ) h jv dz~Q +b (0, z~, 0 T )] in eq. (J2]) stands for the light-like 
Wilson line. Since the light-like Wilson line in QCD involves soft-collinear gluon held 
kD a (x), it can be treated classically as W ta (x) becomes the SU(3) pure gauge when Wil¬ 
son line becomes light-like (see section II for details). Hence for the better understanding 
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of the number operator interpretation of the gluon in eq. (J2J) let us replace the light¬ 
like Wilson line 'PexpfzgT^ 4 ^ fg dz~Q +b (0, z~, Ot)] in eq. (J2J) by its classical counterpart 
Pexp [igT^h dz~A +b (0,z~,0 T )} where A^ a (x) is the classical SU(3) pure gauge back¬ 
ground held. With this one finds from eq. (J2D that F 9 ' l ' a (x)F^ v { 0) in eq. ([2]) contains cubic 
and quartic powers of the (quantum) gluon held Q^ a (x) which implies that the present def¬ 
inition of gluon distribution function at high energy colliders in eq. (J2j) is not consistent 
with the number operator interpretation in quantum mechanics, he., F^ wa (x)F^ u (0 ) is not 
proportional to quadratic power of the (quantum) gluon held Q^ a (x). This implies that the 
present definition of the gluon distribution function at high energy colliders is not correct. 

In this paper we derive correct definition of the gluon distribution function at high energy 
colliders from hrst principles which is consistent with the number operator interpretation of 
the gluon and is gauge invariant and is consistent with the factorization theorem in QCD. 

We find that the correct definition of the gluon distribution function at high energy 
colliders which is consistent with the number operator interpretation of the gluon and is 
gauge invariant and is consistent with the factorization theorem in QCD is given by 

f g /p ( x ) = ^j dy~e~ ixP+ y-x < P\Q;(0,y-,0 T )[Vex P [igT^ b J o V dz~A +b (0, z~ , 0 T )]]Q^(0)|P > 

(4) 


which is valid in covariant gauge, in light-cone gauge, in general axial gauges, in general 
non-covariant gauges and in general Coulomb gauge etc. respectively. Note that the correct 
definition of gluon distribution function at high energy colliders in eq. ([4]) is gauge invariant 
with respect to the gauge transformation 


T“T fl (x) = U(x)T a A a Jx)U~\x) + -[d^U^p-^x), 


W 


along with the homogeneous transformation 


U(x) = e i9Ta “ a{x) (5) 


11 


14] 


T a Q'Vx) = U(x)T a Q a (x)U~\x), 


U(x) = e 


— JgT a u a {x) 


( 6 ) 


where Q^^x) is the hard (quantum) gluon held whose distribution function we measure 
and A tm (x) is the soft-collinear gluon held which is the SU(3) pure gauge background held. 
After renormalization the gluon distribution function is expected to obey a QCD evolution 
equation, like DGLAP equation [ljj], which follows from renormalization group equation. 
We will provide a derivation of eq. (J3J) in this paper. 
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The paper is organized as follows. In section II we describe soft-collinear divergences and 
light-like Wilson line in QCD. In section III we discuss the non-perturbative gluon correlation 
function and proof of factorization theorem in covariant gauge. In sections IV, V, VI and 
VII we discuss the non-perturbative gluon correlation function and proof of factorization 
theorem in general axial gauges, in light-cone gauge, in general non-covariant gauges and in 
general Coulomb gauge respectively. In section VIII we derive the correct definition of the 
gluon distribution function at high energy colliders as given by eq. (j4]) . Section IX contains 
conclusions. 


II. SOFT-COLLINEAR DIVERGENCES AND LIGHT-LIKE WILSON LINE IN 
QCD 


Note that in eq. (jT]) the collincar and soft divergences can occur in the Feynman diagrams 
in partonic level cross sections. Hence it is important to show that these collincar/soft diver¬ 
gences are factorized into the definition of the parton distribution function/fragmentation 
function. This is done by supplying Wilson line in the definition of parton distribution func¬ 
tion/fragmentation function which also makes them gauge invariant. Hence we will briefly 
discuss the issue of gauge invariance and the Wilson line for soft-collinear divergences in 
QCD in this section. 

Before proceeding to the issue of gauge invariance and the Wilson line for soft-collinear 
divergences in QCD let us first discuss the corresponding situation in QED. The gauge 
transformation of the Dirac field of the electron in QED is given by 

if/(x) = e iea ^V(z)- (7) 


Hence we can expect to address the issue of gauge invariance and factorization of soft- 
collinear divergences in QED simultaneously if we can relate the cj(x) to the photon field 
A M (x). This is easily done by using Eikonal Feynman rules in QED for soft-collinear diver¬ 
gences which can be seen as follows. 

The Eikonal propagator times the Eikonal vertex for a photon with four momentum k M 


interacting with an electron moving with four momentum in the limit \k\ 
given by [8, 16-24] 


p ■ k + ie 


\p\ is 


(8) 
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In QED the soft divergences arise only from the emission of a photon for which all compo¬ 
nents of the four-momentum k^ are small (\k\ —> 0) which is evident from eq. ((8]). From eq. 
© one also finds that when 0 < |A;| <<< \p\ and p is parallel to k one may find collinear 
divergences. 

However, since 


p ■ k = pok 0 — p ■ k — 


'i + ^ - cose) 

pZ 


(9) 


one hnds that the collinear divergences does not appear in QED because of the non-vanishing 
mass of the electron, i. e., m ^ 0. From eqs. © and © one hnds that the collinear 
divergences appear only when m — 0 and 6 = 0 where 9 is the angle between p and k. 
Since gluons are massless and the massless gluons interact with each other one hnds that 
the collinear divergences appear in QCD. Since a massless particle is always light-like one 
hnds that the soft-collinear divergences can be described by light-like Wilson line in QCD. 

For light-like electron we find from eq. © 


pH 




= e 


( 10 ) 


p ■ k + ie l ■ k + ie 

where is the light-like four-velocity (|/|=1) of the electron. Note that when we say the 
” light-like electron” we mean the electron that is traveling at its highest speed which is 
arbitrarily close to the speed of light (|/| ~ 1) as it can not travel exactly at speed of light 
because it has finite mass even if the mass of the electron is very small. Hence we find that 
if U l is light-like four velocity then the soft-collinear divergences can be described by the 
Eikonal Feynman rule as given by eq. TO- 
From eq. (ITU]) we find 


d A k l-A(k) 
(27r) 4 l ■ k + ie 


= —ei 


dX 


d A k 

W)‘ 


■ A(k) = ie / d\l ■ A(l\) (11) 


where the photon held H M (x) and its Fourier transform A^{k) are related by 

d A k 


A ‘“M = / 


( 12 ) 


Now consider the corresponding Feynman diagram for the soft-collinear divergences in QED 
due to exchange of two soft-collinear photons of four-momenta and k%. The corresponding 
Eikonal contribution due to two soft-collinear photons exchange is analogously given by 

2 r d A k± d A k 2 l ■ A(k 2 )l ■ A(ki) 

J (27r) 4 (27r) 4 (l ■ {k\ + k 2 ) + ie)(l ■ k\ + ie) 
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(13) 


2 -2 
= e i 


/ dX 2 d\\l ■ A(lX 2 )l ■ A(l\i) 
'o J a 2 


ci r °° r°° 

— / dX 2 dX x l-A{lX 2 )l-AilXJ. 
2! jo jo 


2 ! 

Extending this calculation up to infinite number of soft-collinear photons we find that the 
Eikonal contribution for the soft-collinear divergences due to soft-collinear photons exchange 
with the light-like electron in QED is given by the exponential 

e ieJ o °° 

where is the light-like four velocity of the electron. The Wilson line in QED is given by 


ie f x f dx^Aa(x) 
O J XA ^ 


(15) 


When A^(x) = A 11 (IX) as in eq. (1T4|) then one finds from eq. (fT5j) that the light-like Wilson 
line in QED for soft-collinear divergences is given by 251 


J,e J* dx^A^x) _ —ie dXl-A(x+lX) ie dXl-A(lX) 


(16) 


Note that a light-like electron traveling with light-like four-velocity produces U(l) pure 
gauge potential A^(x) at all the time-space position x M except at the position x perpendicular 
to the direction of motion of the electron (l-x = 0) at the time of closest approach 16.1261127]. 
When A /J '(x) = A ti (Xl ) as in eq. (mi) we find 1 ■ f = AM = A ^ 0 which implies that the 
light-like Wilson line finds the photon field A fl (x ) in eq. CHI) as the U(l) pure gauge. The 
U(l) pure gauge is given by 


A^(x) = d^u(x) 


(17) 


which gives from eq. (fl6l) the light-like Wilson line in QED for soft-collinear divergences 


ieoj(x) —ieu)( 0) ie f* dx^Au(x) —ie f°° dXl-A(x+lX) ie f°° dXl-A(lX) 


(18) 


which depends only on end points 0 and x M but is independent of the path. The path 
independence can also be found from Stokes theorem because for pure gauge 


F^(x) = d fi A v (x) - d v A»(x) = 0 


(19) 


which gives from Stokes theorem 


1 § C dxAA l j(x) __ gie dx^dx 1 'F^(x) _ ^ 


( 20 ) 
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where C is a closed path and S is the surface enclosing C. Now considering two different 
paths L and M with common end points 0 and x M we find 


e ie §c dx>1 A »( x ) — e ie f L dx ,1 A /1 (x)-ie / M dx^A^x) _ ^ 


( 21 ) 


which implies that 


e ie So dx>1A A x ) 


( 22 ) 


depends only on end points 0 and x M but is independent of path which can also be seen from 
eq. (1151) . Hence from eq. (1151) we find that the abelian phase or the gauge link in QED is 
given by 


^ —ie d\l-A(x+lX) _ ^ieuj(x) 


(23) 


From eqs. (JTj) and (1551) one expects that the gauge invariance and factorization of soft- 
collinear divergences in QED can be explained simultaneously. 

One can recall that the gauge invariant greens function in QED 


rx 2 

G(xi, x 2 ) —< ip(x 2 ) x expfie / dx^A^x)] x ip(xi) 

J X\ 


> 


(24) 


in the presence of background field A^(x) was formulated by Schwinger long time ago 28]. 
When this background field A^(x) is replaced by the U(l) pure gauge background field as 
given by eq. (1171) then one finds by using the path integral method of QED that 17 


e ieu(x 2 ) < ^(^ >a e * ew Oi) =< ijj[x 2 ) 'il’(xi) > 

/* OO n QQ 

0 dXl ' A ^+ lx ) < $( X2 ) >A e ie Jo w-Mxi+m ( 25 ) 


which proves the gauge invariance and factorization of soft-collinear divergences in QED 
simultaneously. In eq. (1551) the < ^{x?) > is the full Green’s function in QED and 

< ijj(x 2 ) i)(xi) >a is the corresponding Green’s function in the background field method of 
QED. This path integral technique is also used in []j”J to prove factorization of soft-collinear 
divergences in non-equilibrium QED. 

Hence we find that the gauge invariance and factorization of soft-collinear divergences in 
QED can be studied by using path integral method of QED in the presence of U(l) pure 
gauge background field. Therefore one expects that the gauge invariance and factorization 
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of soft-collinear divergences in QCD can be studied by using path integral method of QCD 
in the presence of SU(3) pure gauge background held. 

Now let us proceed to QCD. The gauge transformation of the quark held in QCD is given 
by 

iP\x) =e i9Tauja{x ^(x). (26) 


Hence one hnds that the issue of gauge invariance and factorization of soft-collinear di¬ 
vergences in QCD can be simultaneously explained if uo a (x ) can be related to the gluon 
held A^ a [x). This is easily done by using Eikonal Feynman rules in QCD for soft-collinear 
divergences which can be seen as follows. 

The Eikonal propagator times the Eikonal vertex for a gluon with four momentum k M 
interacting with a quark moving with four momentum p M in the limit \k\ <<< \p\ is given 


by [8j, |ia-24] 


gT° 




(27) 


p ■ k + ie 

In QCD the soft divergences arise only from the emission of a gluon for which all components 
of the four-momentum k M are small (\k\ —)■ 0) which is evident from eq. (l27lh From eq. (1271) 
one also hnds that when 0 < \k\ «< \p\ and p is parallel to k one may find collinear 
divergences in QCD. 

However, since 


_* _* / vn A 

p- k = poh - p-k = \p\\k\(Jl + -p - cos 8) 


(28) 


one hnds that the collinear divergences does not appear when quark interacts with a collinear 
gluon because of the non-vanishing mass of the quark, i. e., m 7 ^ 0 even if the mass of the 
light quark is very small. From eq. (l28j) one hnds that the collinear divergences appear only 
when m = 0 and 6 = 0 where 6 is the angle between p and k. Since gluons are massless 
and the massless gluons interact with each other one hnds that the collinear divergences 
appear in QCD. Since a massless particle is always light-like one hnds that the soft-collinear 
divergences can be described by light-like Wilson line in QCD. 

For light-like quark we find from eq. f[2T]) 

in 


gT a 


p> 


p ■ k + ie 


— gT a 


l ■ k + ie 


(29) 
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where is the light-like four-velocity (|/|=1) of the quark. Note that when we say the ”light- 
like quark” we mean the quark that is traveling at its highest speed which is arbitrarily close 
to the speed of light (|Z| ~ 1) as it can not travel exactly at speed of light because it has finite 
mass even if the mass of the light quark is very small. On the other hand a massless gluon is 
light-like and it always remains light-like. Hence we find that if V‘ is light-like four velocity 
then the soft-collinear divergences in QCD can be described by the Eikonal Feynman rule 
as given by eq. (l29j) . Note that the Eikonal Feynman rule in eq. (|29l) is also valid if we 
replace the light-like quark by light-like gluon provided we replace T fe “ = —if abc . 

From eq. (1291) we find 


gT a 


d A k l ■ A a {k) 

(27t) 4 l ■ k + ie 


= -gT a i 


d\ 


d A k 


e“ * A ( • A a (k ) = igT a / d\l ■ A a {l\) 


(30) 


where the gluon held A ,ja {x) and its Fourier transform A^ a (k) are related by 

d A k 


A» a (x) = 


(2U 


r ^"“(fc)e 


ik-x 


(31) 


Note that a path ordering in QCD is required which can be seen as follows, see also 10]. 


The Eikonal contribution for the soft-collinear divergences in QCD arising from a single 
soft-collinear gluon exchange in Feynman diagram is given by eq. fl30jh Now consider the 
corresponding Feynman diagram for the soft-collinear divergences in QCD due to exchange 
of two soft-collinear gluons of four-momenta /of and k 2 . The corresponding Eikonal contri¬ 
bution due to two soft-collinear gluons exchange is analogously given by 


2 r d A h d A k 2 T a l • A a (k 2 )T b l • A b (k l ) 

® J (27r) 4 (27 t) 4 (l ■ (ki + k 2 ) + ie)(l ■ k\ + ie) 

poo poo 

= g 2 i 2 / d\ 2 / d\iT a l ■ A a (l\ 2 )T b l ■ A b (l AQ 

Jo J a 2 

^2-2 poo poo 

= / d\ 2 / dXiT a l ■ A a {lX 2 )T b l ■ A b (l AQ (32) 

2! Jo Jo 

where V is the path ordering. Extending this calculation up to infinite number of soft- 
collinear gluons we find that the Eikonal contribution for the soft-collinear divergences due 
to soft-collinear gluons exchange with the light-like quark in QCD is given by the path 
ordered exponential 


V exp 


poo 

[ig dXl ■ A a {lX)T a 


(33) 
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where A is the light-like four velocity of the quark. The Wilson line in QCD is given by 


Ve 


ig dx»A“(x)T“ 


which is the solution of the equation 



d^Six) = igT a A“(x)S(x ) 

with initial condition 


(34) 


(35) 


S(xi) = 1 . 


(36) 


When A fia (x) = AD a (/A) as in eq. (j33j) we find from eq. (T34l) that the light-like Wilson line 


in QCD for soft-collinear divergences is given by 


251 


Ve ig f g x dx»AZ(x)T“ = f Ve -ig /“ dAM«(*+IA)T“l pjg J™dMA*(l\)T° 


(37) 


A light-like quark traveling with light-like four-velocity IA produces SU(3) pure gauge 
potential A^ a (x) at all the time-space position xA except at the position x perpenc 
the direction of motion of the quark (l ■ x — 0) at the time of closest approach 16 
When A^ a {x) = A^ a (\l) as in eq. (l33l) we find I ■ f = AM = A ^ 0 which implies that the 
light-like Wilson line finds the gluon field A^ a (x) in eq. (|33l) as the SU(3) pure gauge. The 
SU(3) pure gauge is given by 

T a A a Ax) = —\dM(x)\ U~\x), U(x) = (38) 

ig 

which gives 


icular to 



U(xf) = Ve“ S ’! ta “ A ^ )T ‘u( Xi ) = e i,TV(,,) 


(39) 


Hence when A^ a {x) = AD a (A/) as in eq. (1331) we find from eqs. (137)1 and (TTHl) that the 
light-like Wilson line in QCD for soft-collinear divergences is given by 


Ve ie fo ^ A »( x ') T ' a = e i 9T a ^{x) e ~igT a ui a ( 0 ) 


jy e -ig f o °° dXl.A a (x+l\)T a 


Ve igf~d\vA“{l\)T a 


which depends only on end points 0 and x^ but is independent of the path. The path 
independence can also be found from the non-abelian Stokes theorem which can be seen as 
follows. The SU(3) pure gauge in eq. (l38|i gives 


F; U [A] = dgAl{x) - d v Afa) + gf abc A b fl (x)A c u (x) = 0. 


(41) 
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Using eq. (1411) in the non-abelian Stokes theorem 30] we find 


Ve ig $ c dx»AXz)T° = Vexp [ lg dx^dx V 

J s 


Ve 


ig r dx ,x A b x (x')T b 


K^)T a 


Ve ig£dx" s A°(x")T° 


= 1 
(42) 


where C is a closed path and S is the surface enclosing C. Now considering two different 
paths L and M with common end points 0 and x 9 we find from eq. 


(43) 


Ve ig§ c dx»A°{ X )T a = p exp [ig J^dx»A“(x)T a - ig J^dx^A a ^x)T a 
= [Ve i9 ^L 


dx^A^^T 11 


V(T i9 i M dxVAb ^ Tb 


which implies that the light-like Wilson line in QCD 

Ve ig J*dx»A“(x)T“ 


(44) 


depends only on the end points 0 and x 9 but is independent of the path which can also be 
seen from eq. (1401) . Hence from eq. (1401) we find that the non-abelian phase or the gauge 
link in QCD is given by 


< 3 )^ _ 'Pe~ i9 fo° d u- Aa ( x + l *) Ta — e igT a ui a ( x ) 


(45) 


Note that from eq. (14T| we find the vanishing physical gauge invariant held strength 
square F ,iva [H] [A] when A^ ia (x) is the SU(3) pure gauge as given by eq. 


. Hence in 

classical mechanics the SU (3) pure gauge potential does not have an effect on color charged 
particle and one expects the effect of exchange of soft-collinear gluons to simply vanish. 
However, in quantum mechanics the situation is a little more complicated, because the 
gauge potential does have an effect on color charged particle even if it is SU(3) pure gauge 
potential and hence one should not expect the effect of exchange of soft-collinear gluons 


to simply vanish 


16]. This can be verified by studying the non-perturbative correlation 


function of the type < 0\ijj(x)'iJj(x')'ijj(x")^(x ,,, )...\0 > in QCD in the presence of SU(3) pure 
gauge background held. 

Under non-abelian gauge transformation given by eq. (J5]) the Wilson line in QCD trans¬ 
forms as 


Ve 


ie f Xf dx»A’°-(x)T a TT/ . 

Jx i M = U(x f ) 


Ve 


ie f*f dx^A“(x)T a 


U 1 (x i ). 


(46) 
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From eqs. (1401) and (1461) we find 

Pe -.,/ 1 "«.(, ttt r = u(x)Ve -isS~MA‘ u(x) = exp [ j 9 rV (x)] 

(47) 

which gives from eq. (|45|) 

&{x) = U(x)$(x), = ^(x)U~ 1 (x). (48) 

In the adjoint representation of SU(3) the corresponding path ordered exponential is 
given by 

M®) = Pexp [-ig f°° dXl ■ A c (x + IX)T^ C ] = e ^ A)c “ c W , (T^ c ) ab = ~if abc . 

Jo 

(49) 


To summarize this, we find that the soft-collinear divergences in the perturbative Feyn¬ 
man diagrams due to soft-collinear gluons interaction with the light-like Wilson line in 
QCD is given by the path ordered exponential in eq. d33j) which is nothing but the non- 
abelian phase or the gauge link in QCD as given by eq. (I45j) where the gluon held A^ a (x) 
is the SU(3) pure gauge, see eqs. (l38]h fl39lh (1401) . This implies that the effect of soft- 
collinear gluons interaction between the partons and the light-like Wilson line in QCD can 
be studied by putting the partons in the SU(3) pure gauge background field. Hence we 
find that the soft-collinear behavior of the non-perturbative correlation function of the type 
< 0\'ip(x) / il>(x')'iij(x")'ip(x" , )...\0 > in QCD due to the presence of light-like Wilson line in 
QCD can be studied by using the path integral method of the QCD in the presence of SU(3) 
pure gauge background held. 

It can be mentioned here that in soft collinear effective theory (SCET) 3jJ it is also 


necessary to use the idea of background fields 
distinct gluon fields jlT 


111] to give well defined meaning to several 


Note that a massive color source traveling at speed much less than speed of light can 
not produce SU(3) pure gauge held hi], 2(], 27]. Hence when one replaces light-like Wilson 
line with massive Wilson line one expects the factorization of soft/infrared divergences to 


break down. This is in conformation with the finding in 


32] which used the diagrammatic 


method of QCD. In case of massive Wilson line in QCD the color transfer occurs and the 
factorization breaks down. Note that in case of massive Wilson line there is no collinear 
divergences which is explained in eq. 
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III. NON-PERTURBATIVE GLUON CORRELATION FUNCTION AND PROOF 
OF FACTORIZATION THEOREM IN COVARIANT GAUGE 


Since gluon distribution function inside the hadron is a non-perturbative quantity in QCD 
it is natural to use path integral method of QCD to study its properties from first principles. 
The generating functional in the path integral method of QCD is given by {ll, 33] 


z[J,v,v\ = 


det( 


S(d^ ia ) 

Suj b 


) i^V^j 2c 


[Q]-hi d ^ a ) 2 +^\ i ^ d ^- m +9T a ^Q a M+J-Q+^+^ r i\ 


(50) 


where is the external source for the quantum gluon held Q^ a (x) and r/j(x) is the 

external source for the Dirac held of the quark and 


FL [<a = d,Qt(x) - dMx) + ar bc Q i Jx)Qi(x), 


pa2 

flU 


[<?] = F^[Q]FL[Q], (51) 


determinant det ^^^T ~ ) can be expressed in terms of 
33] but we will directly work with the Faddeev-Popov 


Note that the Faddeev-Popov (F-P] 
path integral over the ghost helds 
(F-P) determinant det( ^%j^ - ) in this paper. The non-perturbative correlation function of 
the type < O| , 0(a;i)'0(^2)|O > in QCD is given by 17] 

1 


< O|'i/’(Xi)'0(.T 2 ) |0 > = 


m 


'tp(xi)i;(x2) 


Xdet e * f d 4 x[-iF^JQ]--L(d / xQn 2 +^d / x-m+gT^Q^] 

V 50 J b ’ 


(52) 


Similarly the non-perturbative gluon correlation function of the type < 0|Q“((ri)Q])(.T2)|0 > 


in QCD is given by 


17 j 


< 0\Q;( Xl )Q b u (x 2 )\0 >= J[dQ][d$\m q;( Xi )qI(x 2 ) 

Xdet (^3^1) g* / ^xl-^F^lQj-^^Q^f+^dfx-m+gT^Q^} 
V 5oU b ’ 


(53) 


where the suppression of the normalization factor Z\ 0] is understood as it will cancel in the 
final result (see for example, eq. (I93]l ). 

We have seen in section II that the soft-collinear behavior of the non-perturbative (gluon) 
correlation function due to the presence of light-like Wilson line in QCD can be studied by 
using the path integral method of the QCD in the presence of SU(3) pure gauge background 
held. Hence in order to derive correct definition of the gluon distribution function at high 
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energy colliders we use the path integral formulation of the background held method of QCD 
in the presence of SU(3) pure gauge background held as given by eq. 


Background held method of QCD was originally formulated by ’t Hooft 12] and later 


extended by Klueberg-Stern and Zuber 


1- 


14 ] and by Abbott llj. This is an elegant for¬ 


malism which can be useful to construct gauge invariant non-perturbative green’s functions 


in QCD. This formalism is also useful to study 


chromo held 


quark and gluon production from c 
34J] via Schwinger mechanism [35h to compute /3 function in QCD 


perform calculations in lattice gauge theories 


assical 
36], to 


constant in the presence of chromoheld 


38]. 


37] and to study evolution of QCD coupling 


The generating functional in the path integral formulation of the background held method 


of QCD is given by 


11 


m 


Z[A, J, 7 ], rf\ = 


A ,( SGa (Q\ 


J f d i x[-^F a l„[A+Q}-A-(G a (Q)) 2 +^d l ,-m+gT a ^(A+Qr^+J-Q+f)^+H 


(54) 


where the gauge fixing term is given by 


0"i0) = a„Q<“ + gP'-AAr = /),,[4]Q"“ 


(55) 


which depends on the background held A tJ,a (x) and 


f;„\a + Q] = d r \Al + «] - d„[Al + q; ] + gf*°[A b „ + Q‘]K + «]• 


(56) 


We have followed the notations of 


11 


13] and accordingly we have denoted the quantum 


gluon held by Q^ a and the background held by Note that in the absence of the 

external sources the SU(3) pure gauge can be gauged away from the generating functional 
in the background held method of QCD. However, in the presence of the external sources the 
SU(3) pure gauge can not be gauged away from the generating functional in the background 
held method of QCD. 

The non-perturbative correlation function of the type < 0|^(xi)'i/’(x2)|0 >yi in the back¬ 
ground held method of QCD is given by jl3] 

1 


< oiv^iM^lo >a= 


m 




X 


det ( ) e * / d 4 x[-\F a % v [A+Q]-A- (G“ (Q^+^dv-m+gT*^ (A+Q)“] 

V 50J b ’ 


(57) 
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Similarly the non-perturbative gluon correlation function of the type < Q\Q a {x\)Q b v {x 2 )\Q >a 
in the background field method of QCD is given by 17] 


< 0|Q“(xi)Qt(® 2 )|0 >^= j[dQ][d^][d^\ Q a ^x 1 )Q b l/ (x 2 ) 


x 


defr ay / ^xl-^F^lA+Qj-^iG-iQ^+i’li^^-rn+gT-^iA+Q)^] ( 5 g\ 

v 5u b ’ K 


where the suppression of the normalization factor Z[ 0] is understood as it will cancel in the 
final result (see for example, eq. (193]) ). 

The gauge fixing term G a (Q )) 2 in eq. (1541) [where G a (Q ) is given by eq. (155]) ] is 

invariant for gauge transformation of 


&4“ = gf abc A b ^u> c + <9 M cu a , (type I transformation) 
provided one also performs a homogeneous transformation of Q a [ill, 13]: 


(59) 


SQl = gf Q„u e . 


(60) 


The gauge transformation of background field A a as given by eq. (|59l) along with the 
homogeneous transformation of Q l \ in eq. (1601) gives 

+ Q‘,) = gf^K + «>' + d ^ a ( 6 !) 

which leaves — | F a2 fW [A + Q] invariant in eq. (T54|) . 

For fixed A a , i.e., for 


5A “ = 0, (type II transformation) 


the gauge transformation of Q c f ] 



5Q; = gf abc (A\\ + Q > c + d,u /* 


(62) 


(63) 


gives eq. f[6lT) which leaves —\F a2 [A + Q\ invariant in eq. (154|) . Hence whether we use type 
I transformation [eqs. (]59l) and (1601) ] or type II transformation [eqs. (1621) and (J63])] we will 
obtain the same equation (19T1) . 

n 

It is useful to remember that, unlike QED [17| , finding an exact relation between the 
generating functional Z[J, rj, fj) in QCD in eq. (150]) and the generating functional Z[A, J, g, g] 
in the background field method of QCD in eq. (154]) in the presence of SU(3) pure gauge 
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background field is not easy. The main difficulty is due to the gauge fixing terms which are 


different in both the cases. While the Lorentz (covariant) gauge fixing term — ^(d^Q fJ ' a ) 2 in 
eq. CT in QCD is independent of the background field A tm (x), the background field gauge 
fixing term — ^(G a (Q)) 2 in eq. (15T|) in the background field method of QCD depends on 
the background field A^ a (x) where G a (Q) is given by eq. ([55]) 11 13]. Hence in order to 
study non-perturbative correlation function of the gluon in the background field method of 
QCD in the presence of SU(3) pure gauge background field we proceed as follows. 

By changing the integration variable Q —> Q — A in the right hand side of eq. (I58]l we 
find 


< 0|Q£(zi)Q*(z 2 )|0 >a= j[dQ][difi\[diJ>] [Q“(x i) - Afa^Qlfa) - A b „(x 2 )] 

xdet (^1^91) 1 / d4x l-i Fa UQ]-^( G J(Q)) 2 +^ d »- m +9T a rQZm 
v 5oj b ’ 

where from eq. (f55T) we find 


G a f (Q) = d^ a + gf abc A b ^ c - d^ a = D lt {A](r a - d tJ A IM 


and from eq. (1601) [by using eq. (1591) ] we find 


(65) 


*Q a u = -gf u b Q c u + d,io 


( 66 ) 


The eqs. (164|) . (|65|) and (l66]l can also be derived by using type II transformation which can 
be seen as follows. By changing Q —> Q — A in eq. (158]1 we find eq. (T64|) where the gauge 
fixing term from eq. (|55l) becomes eq. (T65|) and eq. (I63|) [by using eq. (162]) ] becomes eq. 
(|66li . Hence we obtain eqs. (l64li . (|65li and (|66li whether we use the type I transformation or 
type II transformation. Hence we find that we will obtain the same eq. (19T1) whether we use 
the type I transformation or type II transformation. 

Changing the integration variable from unprimed variable to primed variable we find from 
eq. (1641) 


< 0|Q“(xi)Q"(x 2 )|0 > a = J [dQ'][d^'][d^'} [Q'^x Q - A^x^lQ*(x 2 ) - A b u (x 2 )} 

' 6 Ga f(Q') ) e i f d 4 x{-±F“2jQ']-±(G?m) 2 +ft[i 7 »dv-m+gT“rQ'“W] 


xdet(- 


5u b 


(67) 


This is because a change of integration variable from unprimed variable to primed variable 
does not change the value of the integration. 
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The equation 


Qt(x) = Qt M + gf^(x)Ql(x) + B^(x) 


( 68 ) 


in eq. (I66|) is valid for infinitesimal transformation (a; << 1) which is obtained from the 
finite equation 

T a Q';{x) = U{x)T a Ql{x)U~\x) + -[d lx U(x)}U~ 1 (x), U{x) = e igTa ^ x \ (69) 

ig 

Simplifying infinite numbers of non-commuting terms we find 


e ~igT b u b (x) jia e igT c u c (x) 


ij 



where 


(70) 


M ab (x) = 


(71) 


Hence by simplifying infinite numbers of non-commuting terms in eq. (169|) [by using eq. 
(ITU and 26]] we find that 


p gM(x) _ 1 

<?'» = [e’ MM UQl(x) + [ U [d^\x)\, M ab {x) = /“*V’(i). (72) 

Under the finite transformation, using eq. fl 72 [1 . we find 

[dQ'\ = [dQ} det [v^] = [dQ\ det[[e 9M{x) }} = [dQ]exp[Tr(ln[e 9M(x) ])] = [dQ\ (73) 
where we have used (for any matrix H ) 


deti7 = exp[Tr(ln/7)]. (74) 

Similarly the fermion fields transform accordingly, see eq. (1261) . Using eqs. (1721) and (126|) 
we find 


W]W) = (#][*/>], ?[ i-fd„ - m + gT*-fQ%W = - m + gT‘rQ‘„]i>, 

KAQ'] = KM- (75) 


Simplifying all the infinite number of non-commuting terms in eq. 
SU(3) pure gauge A 9a (x) is given by [26] 


A^ a (x) = d 9 cj b (x) 


e gM(x) _ ^ 


gM(x) 


ab 


we find that the 


(76) 
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where M a b(x) is given by eq. CD- From eqs. (l76ji and (1721) we find 


Q'» - A‘(x) = [e^UQ^x), M at (x) = f^V(x). 

Using eqs. fj73|l . (175]) and (177]) in eq. (1671) we find 

0|Q“(xi)Qj)(x 2 )|0 > A = j [dQ][di>\[d^] [e 9M{xi) ] ac Ql(x i) [ e 9M{X2) ] bd Q d u (x 2 ) 


(77) 


< 


x 


U f d 4 xl-zF a lAQ}-£(G^Q')f+j]ird„-m+gT“rQ%] 

1 fiw 6 J 


From eq. 


we find 


GJ(Q') = + , l r h 'A l (>‘ ,r - d„A“°. 


(78) 


(79) 


By using eqs. 117211 and -fTTTI'i in eq. lITHll we find 

p gM(x) _ 1 

GJ(C') = dV^UQ^x) + [ U {d^\x)}\ 


+gf abc [d l ‘io e (x) 


e gM{x) _ 

gM ( x ) 


gM(x) 

le 3M{l) UQt(x) + [ 


be 


e gM(x) _ ^ 

gM (x) 


cd [d^ix)}} 


-d^[d 9 uj b (x) 


e gM(x) _ ^ 


gM ( x) 


(80) 


- ab 


which gives 


GJ(Q') = SV^UQ^x)] 

g gM{x) _ ^ 


+gf° bc {dV(x) 
From eq. (T8T1) we find 


gM ( x ) 


p gM(x) i 

]{{e sM,x) UQt(x) + [ U [d,x\x)]]. (81) 


be 


GJ(O') = (f\\e^UQ^ix)} + gr 6c [[9"w'(x)] 
which gives 

G%Q’) = [e^Ud^x) 


gM(x 

e gM(x) _ i 

gM(x) 


\Ue‘ MM UQt(x)} ( 82 ) 


be 


+Q b Jx)&‘{{e° M (% b } + \[dV(x)\ 


g gM(x) _ ^ 

gM (x) 


tif^ti^UQtix)]. (83) 


be 


From 


26J we find 


O^gT^^ij = ig[d^co b (x)} 


e gM(x) _ i 

gM (x) 


<*>]„, M ab (x) = f at v (x) 


ab 


(84) 
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which in the adjoint representation of SU(3) gives (by using T bc = — if abc ) 

' e gM(x) _ i' 


[d^e aM ^} ad = [d^co e (x)j 
Using eq. (J85]) in (l83]i we find 


gM (x) 


gf bac [e M (%d, Mab(x) = f abc u c (x). (85) 


be 


GW) = (x) 


( 86 ) 


which gives 


(G%Q')Y = {d^\x)f 


Since for n x n matrices A and B we have 


(87) 


det(AB) = (detH)(det£>) 


( 88 ) 


we find by using eq. (1861) that 

det [®] = det[ ^l = 

\&ev. mQt{X)) 


det[[e : 


gM(x) I 


5oj b 


= exp[Tr(ln[e 9M(3:) ])] det[ 


5{d^{x)). 

5u b 


= det[ 


Sid.Q^x)). 
5oj b 


Using eqs. (|87li and (|89li in eq. ()78l) we hnd 


< 


x 


0|Q“(xi)Qt(>2)|0 > A = f [dQpVW’] [e 9M{xi) }acQ c ^x i) [e 9M{x2) ] bd Q d u (x 2 ) 

dcU ^QH \ i / d 4 x[-|F“^ [Q] _ l L ( a (iQ ,« ) 2 + ^ [i7 , a(i _ m+9T a / Q^] 

K 8oj b J 


Using the similar technique as above we find 


< 0|>a= 
xdetf ^— e *f d4 A-zF a lAQ}-^(d^) 2 +^d^-m+gT-rQ^ 


<%(*l)<#(* 2 ) 


(89) 


(90) 


(91) 

where 


in the presence of SU(3) pure gauge background field A tJ ' a (x ) as given by eq. 

M ab (x ) is given by eq. (JTTJ). 

Note that eq. (1DTT) is valid whether we use type 1 transformation [eqs. and (1HH ] 

or type II transformation [eqs. (1621) and fj63jl ]. However, since eq. (J5]) is used to study 
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the gauge transformation of the Wilson line in QCD as given by eq. 


, we will use 


type 1 transformation [see eqs. 
transformation give eq. ([5]) and 


and 


0 


1601) in the rest of the paper which for the finite 


13j |. From eq. 


we find that under gauge 


transformation given by eq. © the (quantum) gluon field Q^ a (x) transforms as 


Q'Vx) = [eP M ^UQ l ’lx) 


(92) 


where M ab (x ) is given by eq. (JTTj) . From eqs. (T53l) and (|9T1) we find 




(93) 


From eqs. (1TB . (IBID and (1951) we find 


< 0\Q a (x 1 )Q b u (x 2 )\0 >=< 0\$ IIC (xi)Q c (xi)$ u (x2)Q d t ,(x 2 )\0 > A 


(94) 


which proves factorization of soft-collinear divergences at all order in coupling constant in 
QCD where 

POO 

$ab(x) = Pexp [-ig / dXl • A c (x + l\)T^ c ], (T^ c ) ab = -if abc (95) 

Jo 

is the non-abelian gauge link or non-abelian phase in the adjoint representation of SU(3) 
and is the light-like four velocity. 


IV. NON-PERTURBATIVE GLUON CORRELATION FUNCTION AND PROOF 
OF FACTORIZATION THEOREM IN GENERAL AXIAL GAUGES 


In QCD the generating functional with general axial gauge fixing is given by 


39 


40] 


Z[J,V,v\ = j [dQ] [#][#] 


(96) 


where rf is an arbitrary but constant four vector 


o' 

V 

<M 

Sr 

axial gauge 

v 2 = o, 

light — cone gauge 

to 

V 

o 

temporal gauge. 


(97) 
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Note that unlike covariant gauge in eq. (J50]) there is no Faddeev-Popov (F-P) determinant 


m eq. ijUbj) because the ghost particles decouple in general axial gauges [39|, 40]. The non- 
perturbative gluon correlation function of the type < 0|Q“(xi)Qk(x 2 )|0 > in QCD in general 
axial gauges is given by 


< 0|Q“(zi)Qt(z2)|0 >= 


Q°(xM(x 2 ) 


xe 


i f d 4 xl-iF^[Q]-£(^Q-) 2 +^[i 7 ^-m+gT^Q-M 


(98) 


The generating functional in the background held method of QCD with general axial 


gauge fixing is given by 


40] 


Z[A, J, rpr?] = 

xe* / d4a; [ _ i- Fa ^[ A+< 3]^^D'* ( 3D 2 +V'[^ M9 M- m +9 ra 7' i (^+(3)fi]b+^-Q+^’+bd (gg) 

The non-perturbative gluon correlation function of the type < 0|Q“(xi)Qk(x 2 )|0 >a in 
the background held method of QCD in general axial gauges is given by 

< 0|Q£(xi)<2*(x 2 )|0 > a = j[dQ \[#][#] Q^{xi)Q h v {x 2 ) 

xe i J d 4 x[-\F-l v [A+Q)-^{^Ql) 2 +^d^- m +gT^{A+Q)^}_ Hqq) 


By changing the integration variable Q —> Q — A in the right hand side of eq. (HOOD we 
find 


< 0|Q“Oi)<XO 2 )|0 >a= 


[QuOl) - ^£(xi)][Q*(x 2 ) - ^(x 2 )] 


xe 


i J <i 4 x [ -iF'*^ [ Q ] -i () ,qQ-A ) “ ) 2+vi [ i 7 /‘a M - m + 5 ryQ^] 


( 101 ) 


Changing the integration variable from unprimed variable to primed variable we find from 
eq. (I101D 


< 0|«Ui)«(i 2 )|0 > 4 = / [dQ'][#'][#'] 1<X( X1 ) - A-JxMM - KM] 


xe 


i J d 4 x{-±F a2 M'}-±( V v(Q'-A)Zf+j'{i 7 v d ^- m+gT “ 7 »Q'“}Tp'l 


( 102 ) 


This is because a change of integration variable from unprimed variable to primed variable 
does not change the value of the integration. 

Using eqs. (1751) . ([75]) and (177D in eq. (11021) we find 

< 0|Q“(xi)Q^(x 2 )|0 > A = J [dQ][d$][di/j] [e 9M{xi) ] ac Ql(x 1 ) [e 9M{x2) ] bd Qt(x 2 ) 

xe i f d 4 x[-lF^[Q]-^( v ^([e^ M ^] ab Q b ^x))) 2 +^[i 7 ^-rn+gT a ^Q^] QQ3) 
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which gives 


<0\Q“(x 1 )Q b l/ (x 2 )\0 > A = J[dQ][dip][d'i/}] [e gM{xi) } ac Ql{x i) [e gM(X2) ] bd Q d (x 2 ) 
xe i J d 4 x[-iF“^ [Q] _^ (7?M Qa ) 2 + ^ [i7M ^_ m+gT a rQ a ]V , ] ^ 


(104) 


Using the similar technique as above we hnd 

< 0\[e gMixi \ c Ql l (xi)[e gMix2 \ d Q d u {x 2 )\0 > A = J[dQ][d^>}[dip\ Q a ^(x 1 )Ql(x 2 ) 

xe i f dU[-iF“^[Q ] _^( ?? MQa)2 + ^ [i7 M 5M _ m+gT a 7 MQa] V ,] 

in general axial gauges in the presence of SU(3) pure gauge background field A^^x) as given 
by eq. (I38|) where M ab (x) is given by eq. CD- 
From eqs. (1981) and (11051) we find 

< 0\Ql(xM(x 2 )\0 > = < Q\[e gM ^]acQ C ^( x l)[ e9M ^]bdQt( x 2 )\® >a (106) 

in general axial gauges. From eqs. (1711) . (fF?j) and (11061) we find 

< 0\Q“(xi)Q b v (x 2 )\0 >=< 0\$ac(xi)Ql l (xi)® bd (x 2 )Q d (x 2 )\0 > A (107) 

which proves factorization of soft-collinear divergences at all order in coupling constant in 
QCD in general axial gauges where the non-abelian gauge link or non-abelian phase Q a b(x) 
in the adjoint representation of SU(3) is given by eq. (1951) . 


V. NON-PERTURBATIVE GLUON CORRELATION FUNCTION AND PROOF 
OF FACTORIZATION THEOREM IN LIGHT-CONE GAUGE 


The light-cone gauge corresponds to 


39F41] 


0 ■ Q a = 0, rj 1 = 0 


(108) 


which is already covered by eqs. (PjD and (1971) where the corresponding gauge fixing term 
is given by -^(ri^Q^) 2 . I n the light-cone coordinate system the light-cone gauge 41] 


Q +a = 0 


(109) 


corresponds to 


if = (V + ,V ,V±) = (0,1,0) 


( 110 ) 
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which covers rj ■ Q a = 0 and rj 2 = 0 situation in eq. (11081) . 

Since eq. (11061) is valid for general axial gauges it is also valid for rj 2 = 0 as given by eq. 
(197ft . Hence we find from eq. (llOOp that 

< OlQ^zOQ^IO >=< 0|[e sM ^)] ac g^(x 1 )[e^^] 6d Q^(x 2 )|0 (111) 

in light-cone gauge. From eqs. (171]) . (1791) and (11111) we find 

< 0|Q“(aji)Q*(j? 2 )|0 >=< 0\^ ac (x 1 )Q c fl (x 1 )^ bd (x 2 )Qt(x 2 )\0 > A (112) 

which proves factorization of soft-collinear divergences at all order in coupling constant in 
QCD in light-cone gauge where the non-abelian gauge link or non-abclian phase <f> a b(x) hi 
the adjoint representation of SU(3) is given by eq. (1971) . 


VI. NON-PERTURBATIVE GLUON CORRELATION FUNCTION AND PROOF 
OF FACTORIZATION THEOREM IN GENERAL NON-COVARIANT GAUGES 


In QCD the generating functional with general non-covariant gauge fixing is given by 


42 


43| 




S( n 

det( 1 ’ 


xe 


i f d 4 x[-±F a ? ll/ [Q]-±(2^d /J ,Q2) 2 +i/>[i'y> I d / j-m+gT a 7VQa]i/, + j.Q + . l pr)+r J ip\ 


(113) 


where is an arbitrary but constant four vector. The non-perturbative gluon correlation 
function of the type < Q\Q a ^(xi)Q b v (x 2 )\Q > in QCD in general non-covariant gauges is given 


by 


< 0\Q;(xM(x 2 )\0 >= j[dQ][dip][dip\ Q“(xi)Q^(x 2 ) 


x 


de t ( ^ g) > C-H- 




(114) 


The generating functional in t 
covariant gauge fixing is given by 


re background field method of QCD with general non- 


42 


43] 


,sg a (Q), 


Z[A,J,rj, rj\ = J [dQ] [dip] [dip] det(- , 

Xe »/ d4 ‘ x [-J Fa lv[ A +Q}-^(G a (Q)) 2 +'ip[ i 'y ll d IJ ,-m+g Ta 'yl J -( A +Q)a]i/,+J.Q+fji.’+4’ii] 


(115) 
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where 


rj^rj 1 ' 


<?“(<?) = '^-r (d„Ql + gf^A^Qt) = D„[A}Q‘v 


. . .. .. „. (116) 

r/ r 77 

The non-perturbative gluon correlation function of the type < 0\Q“( X i)Q b (x 2 )\0 in 

the background field method of QCD in general non-covariant gauges is given by 


<o\q;(x 1 )q1(x 2 )\o > a = J[dQ)[dip}[di>\ q;( Xi )qI(x 2 ) 

Suj b 


x 


By changing the integration variable Q —> Q — A in the right hand side of eq. (II 17ft we 
find 


< 0|Q“(ou)Q^(x 2 )|0 >a= 

,mQ). 


[Q a J Xl ) - AV Xl )][Ql(x,) - A b v (x 2 )] 


xdet(—) e*/ d4;r [-2 i;,0 Md<3]-^(6/( , 3)) 2 +^[*7' 1 ^-^+9'r“7' i Qyb] 
5uj b 

where from eq. (Ill6ft we find 


(118) 


rfrf 


Q a f (Q) = (W + gf abc KQt - <V“) - ^ 9f abc (v ■ ^)fa • ^ c ) 

n^n 




(D,[A}Q 


—— r) A a 


(119) 


Changing the integration variable from unprimed variable to primed variable we find from 
eq. (PUP 

< 0|Q“(zi)Q^(x 2 )|0 >a= [[dQ'][d$'][<hl>'] [Q'“(xi) - A a ( Xl )}[Q ,b ( X2 ) - A b u ( X2 )} 


x 


dct( ^^) ) e i / d4;r [-3 Fa MdQ , ]-2k(^(Q , )) 2 +b'[*7' i ^-m+ 9 r“ 7 i‘0^]y']_ ( 120 ) 

5u b 


This is because a change of integration variable from unprimed variable to primed variable 
does not change the value of the integration. 

Using eqs. (1731) . (1751) and (1771) in eq. (112(1 we find 


< 0\Q‘] 1 ( Xl )Q b l/ ( X2 )\0 > A = j[dQ][dip][dip] [e 9M(;Xl) ) ac Q c ^ X i) [ e 9M{X2) ] bd Q d p (x 2 ) 


x 


dctf ^) ) e*/ d4a; [-i J ' 0 Ml c ?]-^^/(Q , )) 2 +V i [n' 1 9M-m+sT“7' i Qyb] 

y 5u b ’ 


( 121 ) 


From eq. (11191) we hnd 


G a f m = -^-[W + gf abc A\ IQ';] - 
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( 122 ) 























By using eqs. (1721) and (1751) in eq. (j!22|) we find 

s p9M(x) _ i 

S/(<?) = + ' 1 ,a ■ 


ry 

+gf°*[d,u*(x) 

I 

2 -^[cW^a; 


gM (x) 


} ab [d u u b (x)}\ 




rj* 


e gM{x) _ ^ 

gM{x) 

e gM(x) _ y 

gM (x) 


e gM{x) _ 2 


+ r gMW ' W 


ab 


which gives 


( 123 ) 


S?(<?) = ^[8J[e»"<*>UQ‘(z)l 


+g/“ fc (a„ W '(i) 

From eq. (1124ft we find 


e gM(x ) _ ^ 


^M(x) 


p gM{x) _ -| 

][[ e’^UQKx) + [ gM(j) U \d^ d (x)}}]. (124) 


be 


which gives 


Gf(Q') = ^IMl^UQU.*)] 


+s/“[[a^ e (x)] 


e gM(x) _ ^ 

gM(x) 


W smx) UQ d M(\ 


be 


(125) 


Using eq. 


G°,m = ^[[e»" (l) ]^<X(i) + <^W4.[[e sM<l) ]a»] 


+[[d„u; e ( a 0] 

in ()126p we find 


e gM(x) _ 1 

gM(x) 


)gJ^\[e^} ed Qi(x)}}. 


- be 


(126) 


which gives 


GKO 1 ) = *B„Qt(x) 


(127) 


,a(n/)) 2 .T/V 


(Gf(Q')r = {^d,Ql(x)Y. 

1 T jZ 

From eqs. (j88l) and (11271) we find 

WJW'), _ T)'‘i7‘'5[[e» M( * ) ]«a p 05(x)] _ , .ijV. 5(^0'(x)), 

6 1 1 “ 11 fot 1 - det ' „2 l C 1“ I 


(128) 


det[[e : 


gM(x) 


■^rrmorM)- 




5u b 


= exp[Tr(l„[e»«W])]det(^ <S(5 '‘ 0 “ W) 


<5a/ 


- a i .i 1 f r f u 5 ( d t^Qu( x )) - 
[ if 5u b ■ 


(129) 
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Using eqs. (1 128p and (1129ft in eq. (11211) we find 


x 


0\Q^(xi)Qu(x2)\0 >a= J [dQ}[d^][d^\ [e 9M(xi) ] ac Q)j(xi) [e 9M{x2) ] bd Q d (x 2 ) 
dct( r ^ ? ^ ) e */(130) 


Using the similar technique as above we find 


< 0|[ e ffM( " l) ] ac Q^(x 1 )[e sM(a:2) ] M Q^(x 2 )|0 > A = J [dQ] [#][#] Q^xQQ^) 

\ jf d4x [-j Fa lAQ}-^i n z3-d fJ ,Qi) 2 +4>[i-yl J 'd fJ ,-m+gT a ^i J 'Q1 l ]ip] fim\ 

in general non-covariant gauges in the presence of SU(3) pure gauge background field A ,ia {x) 
as given by eq. (138ft where M ab (x ) is given by eq. CD- 
From eqs. (II141) and (1131ft we find 


< 0mxM(x 2 )\Q >=< 0|[^ M « I ')]„Q'(i 1 )[e» M <^>] M Q^(x 2 )|0 >., 


(132) 


in general non-covariant gauges. From eqs. (1711) . (I49|) and (1132ft we find 


< 0|Q!(®i)Qt(*2)|0 >=< 0\$ IIC (xi)Q c (xi)$ bd (x 2 )Qi(x 2 )\0 > A 


(133) 


which proves factorization of soft-collinear divergences at all order in coupling constant in 
QCD in general non-covariant gauges where the non-abelian gauge link or non-abelian phase 
$ a ;,(x) in the adjoint representation of SU(3) is given by eq. (j95l) . 


VII. NON-PERTURBATIVE GLUON CORRELATION FUNCTION AND 
PROOF OF FACTORIZATION THEOREM IN GENERAL COULOMB GAUGE 


In QCD the generating functional with general Coulomb gauge fixing is given by 42,(43] 


z{j,n,rj] = / [<ie][#][#] det(dhUAt!Mh) 

xe i f d 4 xl-iF a2 ^[Q]-^(lg^-^l]d M QS) 2 +i,li 7 ^-m+gT a ^Q^+J-Q+^r,+ijr/>] 


(134) 


where 


n‘ u = (1,0, 0,0). 


(135) 
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The non-perturbative gluon correlation function of the type < 0\Q^(xi)Q b (x2)\Q > in QCD 
in general Coulomb gauge is given by 


< 0|<XOi)<^O 2 )|0 >= 


Qu(x M(x 2 ) 


x det( 


TO 


fj,is _ n^n u 


Suj b 




(136) 


The generating functional in the background held method of QCD with general Coulomb 


gauge fixing is given by 42 


Z[A, J, rj, rf\ = 


det( 


£(Q) 

5oj b 


xe 


*/ d 4 a,[-iF“2 1/ [A+Q]-^(e“(Q)) 2 +y[*7 M ^-m+<?T a 7' i (^+<9)M]h+i-Q+^+hd 


(137) 


where 


Mr," 


vrn 


G a (Q ) = \sT - + 9f abc A b ,Ql ) = \sr - '-^-]d,[A]qi 


M n u 


n h n 


n 


n 


(138) 


The non-perturbative gluon correlation function of the type < Q\Q“(xi)Q b (x2)\0 in the 
background held method of QCD in general Coulomb gauge is given by 

< 0|Q“(*TO(* 2 )|0 >a= f [dQ][d#ilj) Q“(xi)Q^(x 2 ) 

x 


detf^Mv f d4 A-j Fa lA A +Q]-^(e a m 2 +^9^-m+gT-r(A+Q) a ^] ( 13 9) 

y Soj b ’ ■ \ ) 


Hence by replacing '-Ar -s. [,g 

we hnd 


n^n 1 ' 


everywhere in the derivations in the previous section 


< 0|Q“(xi)Q^(x 2 )|0 >=< 0\$ ac {x 1 )Q c fl (x 1 )$ bd (x 2 )Qt{ x 2)\0 >a (140) 

which proves factorization of soft-collinear divergences at all order in coupling constant in 
QCD in general Coulomb gauge where the non-abelian gauge link or non-abclian phase 
&ab(x) in the adjoint representation of SU(3) is given by eq. (l95lh 


VIII. CORRECT DEFINITION OF THE GLUON DISTRIBUTION FUNCTION 
AT HIGH ENERGY COLLIDERS 

Under the non-abclian gauge transformation as given by eq. ([5]) we hnd from eqs. (l45lh 
jgp , (l49j) and (1951) that the the non-abclian gauge gauge link or non-abelian phase in QCD 
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in the adjoint representation of SU(3) transforms as 




ab 


= 


0 gM{x) 


<h(x) 


M ab (x ) = f abc cu c (x). 


ab 


(141) 


Hence from eqs. (J92]) and (11411) we find that < 0\^ac( x i)Q^i(xi)^u(x2)Qi(x2)\0 >A in eq. 
(194|) [or in (11071) or in (11121) or in eq. (1 1331) or in eq. (11401) ] is gauge invariant and eq. (194|) [or 
(11071) or (11121) or eq. (I133|) or eq. (11401) ] is consistent with the factorization of soft-collinear 
divergences at all order in coupling constant in QCD. 

Hence from eqs. (194|) [or (11071) or (I112p or (I133p or (11401) ] and (1951) we find that the correct 
definition of the gluon distribution function at high energy colliders which is consistent with 
the number operator interpretation of the gluon and is gauge invariant and is consistent 
with the factorization theorem in QCD is given by 
P+ 


fg/p( x ) = 


x2'K 


dy e 


- -ixP+y 


x 


< 


P I [Ve ~ i9 r dX l^{y+Xl)T^ ]Q a (y) p e -ig / “ rfA l.A°MTW°} QIM (0) |p 


> 


y+=y T =0 

(142) 


which is valid in covariant gauge, in light-cone gauge, in general axial gauges, in general 
non-covariant gauges and in general Coulomb gauge etc. respectively. Eq. (I142p can be 
written as 

x) = ^~ l dy-e~ ixP+y ~x < P\Q a My-,0 T )[Pe i9 C y ’° t) dz - Ac ^ TiA)c ]Q^(0)\P > 
XZ7T J ^ 

(143) 

which reproduces eq. (j4j). This completes the derivation of the correct definition of the 
gluon distribution function at high energy colliders from first principles. 



IX. CONCLUSIONS 

Unlike QED, since F^(x)F^ a (0) in QCD contains cubic and quartic powers of the gluon 
held the present definition of the gluon distribution function at high energy colliders is not 
consistent with the number operator interpretation of the gluon. In this paper we have 
derived the correct definition of the gluon distribution function at high energy colliders 
from Erst principles which is consistent with the number operator interpretation of the 
gluon and is gauge invariant and is consistent with the factorization theorem in QCD. 
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After renormalization the gluon distribution function is expected to obey a QCD evolution 
equation, like DGLAP equation [ijj], which follows from renormalization group equation. 
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